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We investigate the cosmological evolution of gravitational waves in Friedman-Robertson-Walker 
brane world models embedded in a five dimensional anti de-Sitter spacetime. To predict the spec- 
trum of stochastic gravitational background at present generated in the inflationary phase, we 
numerically calculate the evolution of gravitational waves according to the method developed in our 
previous paper [K. Ichiki and K. Nakamura, arXiv:hep-th/0310282 . The resulting logarithmic en- 
ergy spectrum is significantly altered from that of the standard four dimensional infiationary models 



in the frequency range / > /AdS 
de-Sitter spacetime, I ~ 1mm. 



10 Hz, which corresponds to the curvature scale we set of anti 



In recent years much attention has been paid to inter- 
esting models of spacetimes in which fields of the stan- 
dard model are confined to a "brane", while only the 
gravitational field propagates in all the dimensions called 
the "bulk" . Since the proposal of a brane world model of 
our spacetime by Randall and Sundrum 'l'| (RS II model) , 
the phenomenology of brane world cosmological models 
has been the subject of intensive investigations 2]. In 
their models our four dimensional spacetime is embedded 
in the five dimensional anti de-Sitter (AdS5) spacetime 
with curvature scale I, and the four-dimensional gravity 
is effectively recovered in larger scales than I. The cur- 
rent direct bound of the experimental gravitational probe 
allows the scale of extra dimensions or curvature scale / 
to be as large as a (sub)milfimeter |3|. 

In addition to the direct experimental probe of the 
Newton potential, if our universe is really the RS II type 
brane universe, there is a possibility that the evolution 
history of the universe is significantly altered when the 
Hubble scale H~^ is smaller than I in the early universe. 
Actually, some authors have been discussed new physics 
in the epoch of the early universe where Hubble scale is 
smaller than I and the possible sources of gravitational 
wavesjj,l3- In particular, it is pointed out that these new 
physics in the early universe will appear in the stochas- 
tic gravitational wave background with the characteristic 
frequency around 10~^ Hz|j|, if the scale of / is order of 
1mm. 

Thus, the stochastic gravitational waves will be a 
promising candidate which provides direct and deep 
probes of the early universe (for example, see [a])- One 
might see in principle earlier epochs of the universe than 
the photon last scattering surface by gravitational waves. 



In order to give theoretical predictions of the stochastic 
gravitational waves, we have recently proposed the single 
null coordinate system to solve the gravitational waves 
in the expanding brane universe[3|, based on the idea in 
Ref. |8| . We have concentrated on the brane world model 
of the fiat Friedmann- Robertson- Walker (FRW) universe 
without "dark radiation" following discussions in Ref.[9| 
and the fact that the fiat FRW universe is supported 
by recent precise measurements of the cosmic microwave 
background (CMB) anisotropics [l3- Though some au- 
thors have investigated based on the Gaussian normal 
coordinate system in the neighborhood of the brane |llj| , 
this coordinate system includes a coordinate singularity 
in the bulk, where the metric function vanishes [l3| , and 
the treatment of wave equation near the singularity be- 
comes delicate. However, these delicate problems are re- 
moved by the introduction of the coordinate system pro- 
posed in Ref.Q. 

In this letter, we investigate how the presence of extra 
dimension of RS type leaves an imprint on the stochastic 
gravitational wave background at present. We solve the 
evolution of gravitational waves in the brane world model 
and predict the present energy spectrum, following the 
previous proposal|7|- In terms of the single null coordi- 
nate system proposed in Q, the metric on the AdS5 is 
given by 



ds' 



F{T,w)du''-2F{T,w)dudw + r^{T,w)dEo . (1) 



Here the metric functions _F(r, w) and r{T,w) are given 

by 



F{t{u,w),w) = 



ryp 



r(T{u,w),w) ^ -—- X 



^/ryW+a^ + d 

(l-A)e^/' + (l + ^)e-^/' 



(2) 
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respectively, where A = y/l+WH^. Here a{T) is scale 
factor of the FRW brane whose time evolution is given 



by the generalized Friedmann equation [3] 

SttGn A4 



^'^If/^^^ + T + ft"^^ ,3, 



In Eq.Q, a — da/dr, t is the comoving time on 
the FRW brane, p is the energy density on the brane, 
A4 = k|A^/12 + 3A5/4 is the cosmological constant in- 
duced on the brane, Gn = K|A/487r and K5 are the 
four-dimensional and five-dimensional gravitational con- 
stants, respectively. We also assume Z2 symmetry across 
the brane following to the original RS II model. Hereafter 
we set I = 1 mm for the purpose of numerical calcula- 
tions. 

Apart from the polarization tensor, the equation for 
gravitational waves in the bulk is simply given by that 
for the five dimensional massless scalar field, O^h — 0. 
As derived in Ref.[7|, this equation is equivalent to: 



duh 



-3/2 




( dq\ k^ 

dwhiJ~g=]~ g^h ] dw + Ciu) 



(4) 



where — fc^ is the eigen value of the Laplacian of Sq, and 
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U^ ' 2 dw 


2 du 



r.3/2 



(5) 
(6) 



C{u) is determined by the boundary conditions of the 
Neumann type on the brane at each time step 



C{u) = a^/'^F{T,w) X ckjh{u,w)\w=o 



(7) 



Once given an initial spectrum on an initial null hyper- 
surface in the bulk, which is a m = const hypersurface in 
the bulk, Eqs.(|nil-|ZIl are enough to predict the evolution 
of gravitational waves. 

As an initial condition of numerical examples, we fol- 
low the inflationary scenario and we start numerical cal- 
culation by setting h = const, on a initial null hypersur- 
face. This initial value is a solution to the wave equa- 
tion in the long wavelength limit k/aH -^ [l3|- We 
should note, however, that, for fc ^ mode, h = const. 
is not an exact solution to the wave equation but an ap- 
proximate one in the situation k/aH <C 1. Further, as 
the universe expands and decelerates, the approximation 
k/aH <C 1 becomes violated. Conversely, this also im- 
plies that h = const, is an appropriate initial condition 
in the long wavelength approximation, if we start the cal- 
culation from the past much earlier time than that of the 
horizon crossing of each mode on the brane. For this rea- 
son, we start the calculation when the wave is far enough 
outside the horizon with the initial condition h = const. 
to keep the validity of the long wavelength approximation 
k/aH -^ 0. 
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FIG. 1: The time evolution of gravitational waves with co- 
moving frequency, f — 5 x 10^^ [Hz] (upper panel) and 
f — 5 X 10"* [Hz] (lower panel). The time evolution of grav- 
itational wave with / = 5 x 10"^ [Hz] is almost "four dimen- 
sional" one, in which amplitude of gravitational waves scales 
as /i4D oc a~^. On the other hand, the amplitude of gravi- 
tational wave with / = 5 x 10~* [Hz] in brane world model 

decreases faster than a^^ when Aphy ~ /. After that, its 
evolution asymptotically becomes similar to that of standard 



model when A 



phy 



/. 



We show the examples of the evolution of gravita- 
tional waves /i5d in Fig^ which are obtained by solving 
Eqs.Q-Qj numerically. In these figures we also show 
"four dimensional gravitational waves" h^d for compari- 
son, which is a solution of the four dimensional equation 
/i4d + 3-ff/i4d+ (fc^/a^j /i4d = and simply evolves as a~^. 

The feature of the evolutions of gravitational waves, 
which one can easily see from Fig^ is that the ampli- 
tudes of them are constant in time at first, and then 
start damping with oscillation. In the standard four- 
dimensional cosmology, it is well known that the ampli- 
tude of gravitational waves is "frozen" when the modes 
are outside the horizon {k/aH <C 1), and then starts to 
oscillate after the modes enter the horizon {k/aH ^ 1). 
As shown in Fig^ it turns out that this behavior is 
also seen in the evolution of gravitational waves in brane 
world model considered here. 

The other important feature is that the amplitude of 
five-dimensional gravitational waves /isd decreases faster 
than that of standard four-dimensional gravitational 
waves /i4d in higher frequency modes (the lower panel in 
FigQ). This is an essential feature of the five-dimensional 
evolution of gravitational waves. As we mentioned above, 
the gravitational waves with wavelength smaller than the 
curvature scale of AdSs (Aphy = 2'Ka/k < I) would feel 



the universe to be five dimensional, and thus the evo- 
lution of gravitational waves with such wavelength can 
be significantly altered. Even though the decomposition 
of independent Kaluza-Klein (KK) modes is quite non- 
trivial in FRW brane universe, the strongly decelerating 
brane motion and Doppler effects by moving brane cause 
the mode mixing of gravitational waves, and do generate 
the effective "KK modes" which propagate into the bulk 
[ll|. This is the reason for the damping in amplitude 
of fivc-dimcnsional gravitational waves relative to that of 
four-dimensional gravitational waves. 

In Fig[21 we show both the important criterions dis- 
cussed above. We depict the scale factors when each 
Fourier mode with comoving frequency /(= fc/27r) enters 
the Hubble horizon (Aphy = ^1 H; hue), and when their 
wavelength becomes longer than the curvature scale of 
AdSs (Aphys — I] dashed hue) in Fig[21 A break around 
/ « 5 X 10^"* Hz in the line Aphy = \/ H corresponds to 
the frequency of the mode which crosses the Hubble hori- 
zon at the transition from p^-term dominate era to the 
standard radiation dominate era (see, Eq. 0). These 
two lines cross each other (Aphy = ^/ H — I) at the co- 
moving frequency /AdS ~ 10^** Hz. Figure |2] is helpful 
to understand that gravitational waves with comoving 
frequency larger than /acIS Hz are affected by the exis- 
tence of the fifth dimension, and the gravitational waves 
change their behavior at / = /Ads 

To explain this, let us first consider the lower frequency 
modes, / < /Ads- Figure|21shows that the physical wave- 
length Aphy is larger than I first, but the mode is still far 
outside the horizon. After that, as the universe expands, 
the mode crosses the horizon and starts to oscillate. Since 
the physical wavelength is long enough relative to the cur- 
vature length of AdSs at that time, the evolution of gravi- 
tational waves is effectively four-dimensional one (see the 
upper panel in Fig^l. 

On the other hand, the gravitational waves with co- 
moving frequency greater than /AdS enter the Hubble 
horizon, first. Then, they start to oscillate and their 
amplitude decreases faster than a^^ due to the "KK 
mode" mixing discussed above. After that, as the uni- 
verse expands the wavelength of gravitational waves fi- 
nally becomes longer than the curvature scale I of AdSs , 
and the time evolution of gravitational waves asymptoti- 
cally mimics that of the standard four-dimensional FRW 
model, i.e., the amplitude evolves as h^d oc a~^ . The evo- 
lution of gravitational waves in the lower panel of Fig^ 
does show this behavior. 

In order to predict the spectrum of stochastic gravita- 
tional waves at present, we have to know how much am- 
plitude of gravitational waves in five dimensional model 
has escaped from the brane compared to that in four 
dimensional model. To accomplish this, we draw two en- 
velope curves of the evolution of gravitational waves for 
/i5d and /i4d, and numerically observe the ratio h^d/h^d 
after the evolution of gravitational waves in five dimen- 
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FIG. 2: Cosmic scale factors at the brane as a function of co- 
moving frequency when each Fourier mode re-enter the Hub- 
ble horizon (Aphy = ^/H; line) and when the wavelength be- 
comes longer than the AdS5 curvature scale (Aphy — l\ dashed 
line) . 



sional model exhibits asymptotic behavior h ex a~^ . The 
result is depicted in Figj^l Although there seem to exist 
some small scatters, we find that the ratio can be well fit- 
ted by a simple power law above some critical frequency: 



^5d 
/l4d 



\ /AdS J 



for / < /AdS 
for / > /AdS 



(8) 



where we can obtain from FigOlthat /AdS ~ 10~* Hz and 
a is roughly about 0.9. 

The spectrum of stochastic gravitational waves can be 
obtained by multiplying the expected amplitude of four 
dimensional gravitational waves today, ft.4d(''o), by the 
factor (/i5d//i4d) in Eq.Q. Present spectrum of four- 
dimensional gravitational waves /i4d(''o) can be obtained 
as follows. In high energy era when p^ term dominates, 
Hubble parameter scales a.s H oc pr (x a^^ (see Eq.Q) 
leading to l/aH ex a^. Using the fact that the amplitude 
of gravitational waves is almost constant until horizon 
crossing k/aH ~ 1 and decreases as /i4d oc a^^ after that, 
we obtain /i4<i ex fc~^" in high energy era. The same 
arguments lead that /i4d oc k~^ in standard radiation 
dominated era, and /i4d oc k^^ in matter dominated era 
[a|. By combining these considerations, we can predict 
the final spectrum of the stochastic gravitational waves. 

In Fig0|we show the logarithmic energy spectrum of 
gravitational waves, QGw{f)hl, where Ogw(/) oc k'^h'^ 
is the energy density of gravitational waves per logarith- 
mic frequency interval in critical density units, and ho is 
Hubble parameter today in units of 100 km/s/Mpc. We 
find that the spectrum becomes red, ficw^o °^ /~°'*^j 
in the frequency range / > /AdS, contrary to the expec- 
tation that the spectrum would have a spike due to the 
non-standard (strongly decelerating) expansion law |l4| . 

Finally we briefly mention the observational aspects. 
In standard models with CMB quadrupole anisotropics 
measured by COBE, the energy density of gravitational 
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FIG. 3: The amplitude ratio of gravitational waves (hsd) to 
the standard four dimensional amplitude (/i4d oc a~^) deter- 
mined in the asymptotic era (Aphy ^ I)- Damping due to 
the momentum into the extra dimension becomes larger with 
larger comoving frequency. 






[ CMB bound 



LISA 



DECIGO/BBO 



ultimate DECIGO 




f[Hz] 

FIG. 4: The spectrum of stochastic gravitational waves in 
brane world model with I = 1 mm (line) at present. The 
case in which we only include the effect of non-standard ex- 
pansion law is shown by dashed line for comparison. Overall 
amplitude is normalized by CMB bound around the frequency 
10"^* Hz. The noise curves of LISA, practical DECIGO/BBO 
and quantum limited ultimate DECIGO for integration times 
~ 3 yr are also shown in the figure [VA HM llfll l2(i| . 



waves is constrained to ficw ^ 10^ at / ^^ 0.1 Hz in 
standard cosmology. If future CMB experiments could 
directly measure tensor amplitude and its spectral tilt 
(nr) through temperature and polarization anisotropics 
[l5J , one can in principle extrapolate and expect the am- 
plitude of gravitational waves to much higher frequency 
range, even up to 10^^ Hz [l^. On the other hand, in 
the brane world models considered in this paper, it is 
derived from our results that the spectral energy density 
will be much smaller than expected above some critical 
frequency /acIS ~ 10^^ Hz. Unfortunately, stochastic 
radiation from white dwarf binaries prevents us from de- 



tecting inflationary primordial gravitational waves in the 
frequency range 10~^ < / < 10~^ Hz. However, the 
sky would be "transparent" to a primordial signal above 
the frequency / ~ 0.1 Hz by one year observation, if 
the neutron star - neutron star merger rate is around 
R ^ 10^^/yr/galaxy [l3. In this bands, the decihertz 
(10^^ — 10 Hz) interferometer gravitational wave detec- 
tors, such as DECIGO [13] and BBO H^^ have been 
proposed aiming to detect the primordial gravitational 
waves as one of their important goals. Therefore, the 
signal of such brane world cosmological models might be 
detected as a "lack" of gravitational waves in such higher 
frequencies. 
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